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ABSTRACT 

A  knowledge  of  the  interference  between  direct,  surface  and  bottom  refelected 
rays  emanating  from  a  moving  underwater  source  can  be  used  to  determine  the 
source-detector  relative  positions.  Using  simple  ray  tracing,  the  frequency-time 
interference  patterns  for  a  broadband  frequency  source  arising  when  the 
propagating  medium  has  a  piecewise  linear  soimd  speed  profile  are  considered, 
and  are  compared  with  those  for  a  constant  sound  speed.  It  is  found  that  even 
small  deviations  from  a  constant  soimd  speed  can  lead  to  a  noticeably  different 
interference  fringe  placement  and  shape,  and  consequently  to  the  inferred 
source-detector  geometry.  Further,  comparisons  of  model  calculations  to 
experimentally  observed  fringe  systems  demonstrate  the  inadequacy  of  the 
constant  sound  speed  model,  whereas  the  use  of  the  simple  segmented  profile  is 
capable  of  replicating  the  observed  patterns. 
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The  Effect  of  Ray  Curvature  on  Lloyd's  Mirror 
Fringe  Estimates  and  Source  Localisation 


EXECUTIVE  SUMMARY 

This  report  considers  the  effect  a  realistic  sound  speed  profile,  or  variation  of  the 
sound  speed  with  water  depth,  can  have  on  the  passive  detection  and  range 
estimation  of  a  source.  Because  of  reflections  from  the  sea  surface  or  floor,  the 
sound  measured  at  a  hydrophone  can  be  the  result  of  wave  interference  between 
the  directly  transmitted  and  a  reflected  sound  wave,  each  of  which  is  emitted 
from  the  source  (the  Lloyd's  mirror  effect).  This  interference  will  cause  the 
soimd  intensity  measured  at  the  hydrophone  to  be  modified  depending  on  the 
source  position  and  the  frequency  of  the  sound.  Hence  for  a  moving  source,  the 
measured  soimd  intensity  will  fluctuate  with  time  and  the  frequency.  Further, 
for  fixed  source  and  hydrophone  positions  the  variation  will  be  dependent  on 
the  way  in  which  the  sound  propagates  through  the  water,  or  equivalently,  the 
way  the  speed  of  soimd  varies  with  depth  (and  location).  Depending  on  the 
structure  of  the  soimd  speed  profile,  the  variation  may  be  important  at  a  given 
distance,  whilst  any  effects  arising  from  the  structure  of  the  sound  speed  profile 
will  become  more  evident  with  increasing  separation  of  the  source  and 
hydrophone.  The  information  inherent  in  these  sound  intensity  variations 
allows  an  accurate  estimation  of  the  source  range,  or  position  relative  to  the 
hydrophone,  but  is  based  on  an  assumed  model  of  sound  propagation.  It  is 
important  then  to  consider  whether  the  simple  model  of  a  constant  sound  speed 
is  adequate  when  extracting  this  information. 

Predictions  based  on  a  variable,  but  still  simple,  soimd  speed  profile  are 
presented.  Firstly,  a  phenomenological  appraisal  of  the  differences  that  can  arise 
by  the  choice  of  either  the  constant  or  the  simple  variable  sound  speed  models  is 
made.  This  is  based  on  several  realistic  sound  speed  profiles  appropriate  for  the 
Australian  costal  region.  Significant  effects  can  be  obtained,  depending  on  the 
size  of  the  differences  between  these  h\m  types  of  profile,  and  are  found  to  be 
present  for  even  small  ranges  (==0.5-1  km).  Thus,  the  type  of  assumed  model 
can  realistically  be  expected  to  affect  the  extracted  source  positional  information 
even  for  passive  detection,  which  is  generally  only  of  short  range. 

We  then  compare  some  true  data  and  the  predictions  based  on  each  of  the  two 
sound  speed  profile  models.  These  data  show  a  strong  Lloyd's  mirror  effect. 
The  comparisons  indicate  that  a  variable  sound  speed  profile  is  necessary  in 
order  that  a  suitable  description  of  some  of  the  data  is  obtained;  the  constant 
soimd  speed  profile  is  not  capable  of  correctly  reproducing  these  data.  We  also 
find  that  even  when  either  model  provides  an  adequate  description,  the 
estimated  source  positions  differ  by  a  noticeable  amount,  contrary  to  the 
expectation  that  the  estimated  source  positions  should  be  similar. 


In  conclusion,  the  Lloyd's  mirror  effect  can  be  used  to  provide  an  acciurate 
estimate  of  the  source  posihon  relative  to  the  hydrophone,  but  this  estimate  is 
dependent  on  the  assirmed  model  of  sound  propagation.  A  variable  soimd 
speed  profile  can  make  a  significant  difference,  and  is  necessary  in  describing 
some  of  the  experimental  data.  Therefore,  it  is  important  that  the  suitability  of 
the  underlying  model  of  propagation  be  ascertained,  and  in  particular  the 
suitability  of  the  constant  sound  speed  model,  before  the  deduced  source 
position  is  accepted. 
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1.  Introduction 


In  this  brief  report  we  consider  the  effect  of  applying  a  realistic  sound  speed  profile  to 
the  determination  of  Lloyd's  Mirror  interference  fringes,  and  the  subsequent 
localisation  of  a  source  position. 

Given  a  broadband  acoustic  source  located  at  some  range  and  depth  from  a  sensor,  it 
is  possible  to  determine  the  (relative)  location  of  this  soiuce  point  from  the  interference 
between  a  direct  and  a  singly  reflected  beam  (Lloyd's  Mirror  effect)  as  observed  in  the 
received  signal  frequency  spectrum;  the  frequency  spacing  of  interference  minima  or 
maxima  being  dependent  on  the  source /receiver  geometry.  Thus,  the  Lloyd's  Mirror 
(LM)  interference  will  modulate  the  frequency  spectrum.  Moreover,  if  the  source  is 
moving,  then  the  time-frequency  spectral  plot  will  show  these  extrema  as  lines  of 
constant  phase  difference.  These  will  overlap  any  narrowband  constant  frequency 
emissions  from  the  source  and  thus  will  give  a  temporal  modulation  of  the  received 
narrowband  intensity.  Even  in  the  absence  of  broadband  emissions,  the  Lloyd's 
Mirror  effect  can  be  observed  through  this  temporal  modulation  of  the  narrowband 
intensity.  Hence,  it  can  offer  an  important  method  of  determining  the  source  position 
directly  from  the  measured  frequency  spectrum  of  a  single  receiver 

A  common  approach  to  determining  this  position  (relative  to  the  receiver)  is  to 
assiune  that  the  water  column  is  an  isotropic,  homogeneous  medium  with  a  constant 
sound  speed.  In  this  case  a  simple  analytical  expression  can  be  derived  for  the 
frequencies  corresponding  to  interference  extrema.  Whilst  in  some  instances  a 
constant  speed  will  be  appropriate,  more  generally  this  is  not  a  suitable  description  of 
the  soimd  speed  profile.  Here,  the  degree  to  which  a  variable  sormd  speed  profile  can 
affect  the  received  interference  pattern  is  considered. 

We  are  motivated  in  doing  this  by  the  need  for  accurate  estimations  of  the  source 
position  to  be  able  to  determine  absolute  source  sound  pressure  levels  of,  in  particular, 
narrowband  emissions,  and  in  the  case  of  a  moving  source,  temporal  modulations  to 
the  narrow  band  intensity  arising  from  the  LM  interference  effects. 

This  report  is  divided  as  follows.  First,  a  description  is  given  of  the  model 
underlying  the  calculations.  Then  follows  an  overview  of  the  effect  of  introducing  a 
non-constant  sound  speed  profile,  and  a  discussion  of  particular  case  studies.  Finally, 
conclusions  are  presented. 


2.  Model  Description 


For  the  sake  of  simplicity,  and  with  little  loss  of  generality,  consider  a  source  at  a  fixed 
depth,  moving  with  a  constant  velocity,  and  having  some  point  of  closest  approach 
(CPA)  to  a  receiver.  With  a  constant  sound  speed  profile,  rays  propagate  in  straight 
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lines,  and  an  analytical  solution  is  available  for  the  positions  of  the  LM  interference 
fringes  seen  in  a  time-frequency  spectrum.  We  wish  to  maintain  the  advantages  of  an 
analytical  solution,  and  thus  consider  a  piecewise  linear  and  continuous  profile.  The 
range  independence  and  isotropic  nature  of  the  medium  remain  unaltered.  A  general 
analytical  solution  is  not  available  in  this  case,  however  with  a  constant  soimd  speed 
gradient  rays  propagate  along  circular  arcs^l  Thus,  an  analytical  solution  for  a  ray 
trajectory  is  available  for  each  profile  segment,  and  raytracing  can  be  done  to 
determine  numerically  the  interference  pattern  seen  in  a  time-frequency  spectrum. 

Consider  a  single  layer  of  the  model  profile  (figure  1).  A  ray  emitted  horizontally 
will  have  a  radius  of  ciuvature  Pg  determined  by  the  sound  speed  gradient  7^^, 


Po(^m)  = 


CsizJ 

Y 


where  is  the  ray  starting  depth  relative  to  the  layer  starting  depth  (zg)  and 
Cliz  =  Zq)  +  '^ ■  is  the  speed  of  sound  at  this  depth.  For  any  other  ray 
emitted  at  some  angle  0;^,  to  the  horizontal,  the  radius  of  curvature  p  is  greater,  with 
the  centre  of  curvature  lying  at  a  depth  pg  below  the  layer  depth  (figvue  1).  That  is 

Po^ 

cos0„ ■ 


Simple  geometry  then  defines  the  quantities  of  interest.  For  a  ray  starting  at  a  depth 
z^  at  angle  9^  and  terminating  at  depth  at  angle  0;^, ,  the  change  in  range  5R  is 

5i?=p-(sine^,  -sin0^). 


the  change  in  depth  is 


5Z)  =  p-(cos0^-cos0^^) 

=  Po-p-cos0,,, 

and  the  increment  in  path  length  is 

6L  =  p-(0^,-0J. 

The  sign  convention  in  use  is  that  downward  rays  have  positive  angles,  and  > 

Zg  >  0.  To  determine  the  interference  effects  between  two  curved  rays,  we  need  to 
consider  the  transit  time  difference  of  the  two  ray  paths.  The  transit  time  along  a  path 
S  is 
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Z  =  0 


/  /  I 


Figure  1:  Single  layer  ray  geometry.  The  ray  with  radius  of  curvature  p  starts  at  depth  Z,- 
below  the  top  of  the  layer  and  terminates  at  depth  Zj . 


Using 

and 

we  obtain 


z  =  Py  -p-cosG 

=  Q(^,n)+Y-Po-Y-p-cose 
=  -7-p-cos6 

87  = 

Y  ,/  cos  9 

y  ui 

=  —  j  atanh(tan(— ))  -  atanh(tan(^^)) 
Y  I  2  2  . 


These  form  the  basic  equations  with  which  a  ray  can  be  traced  from  a  source  to  a 
receiver.  In  the  calculations  the  interest  is  in  the  interference  effects  arising  from  a 
direct  and  a  singly  reflected  ray.  Hence,  the  aim  is  to  determine  the  transit  times  of 
these  two  rays.  For  the  case  of  broadband  emissions,  the  transit  time  difference  is 
used  in  conjrmction  with  a  desired  phase  difference  to  define  a  frequency  (eg.  the 
frequency  at  which  the  interference  fringe  of  a  given  order  occurs),  whilst  for 
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emissions,  the  frequency  defines  a  phase  difference,  and  thus  for  a  moving  source  a 
temporal  intensity  envelope. 

For  a  surface  reflected  ray  a  n  phase  shift  is  introduced,  and  the  frequencies  f^,  for 
order  n  corresponding  to  interference  extrema  are  given  by 


fr  =  (n  +  l)- 
fr=(^+j)- 


^^irect  ^^eflect 


for  n  =  0, 1, .  •  • 


For  bottom  reflection,  the  extrema  frequencies  are  exchanged.  Each  order  n 
corresponds  to  a  particular  fringe  in  the  time-frequency  spectrum. 

For  later  identification  purposes,  CR  and  SR  will  refer  to  the  curved  and  straight  line 
ray  propagation  models  respectively. 


3.  Lloyd's  Mirror  Calculations 


Two  approaches  to  assessing  the  interference  effects  arising  from  the  CR  model  profile 
have  been  taken.  Firstly,  a  set  of  realistic  profiles  is  considered,  with  the  intention  of 
drawing  some  more  general  conclusions  concerning  the  differences  between  the  SR 
and  CR  models.  Specific  profiles,  together  with  data,  are  then  considered,  and  serve  to 
highlight  differences  between  the  two  models.  To  avoid  complications  that  may  arise 
with  the  bottom  layer,  only  surface  reflections  are  considered.  This  is  a  reasonable 
assumption,  since  for  many  situations  of  interest  the  water  depth  is  large  and  bottom 
bounce  interference  effects  will  be  diminished  via  signal  attenuation  and  will 
correspond  to  very  low  frequency  fringes. 


3.1  Overview  of  CR  Model  Results 


A  set  of  twelve  sound  speed  profiles  has  been  selected,  covering  six  locations  about  the 
Australian  coast,  and  two  seasons  of  the  year.  Each  was  obtained  from  the  TESS^^ 
database,  wherein  sound  speed  profiles  are  gridded  at  1°  intervals.  Each  databased 
profile  is  selected  from  a  set  of  profiles  measured  from  within  the  appropriate  grid 
region,  and  is  assessed  as  best  representing  the  profile  features  from  this  set.  Thus,  the 
TESS  sound  speed  profiles  can  be  expected  to  be  realistic  whilst  giving  the  general 
sound  speed  variations  for  a  particular  region.  The  two  seasons  are  represented  by  the 
January  and  July  data.  The  collection  of  twelve  profiles  represents  the  variation  that 
can  be  expected  in  sound  speed  profiles  in  the  Australian  region.  A  summary  of  these 
is  given  in  Table  1  and  Figure  3. 
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Table  la:  Sound  speeds  and  RMS  sound  speed  deviations  relative  to  the  surface  sound  speed 
as  used  in  the  calculations  described  in  section  3.1.  Data  is  for  the  month  of  January. 


Tasman  Sea 

off  Brisbane 

Coral  Sea 

Exmouth  Gulf 

off  Perth 

Great  Aust.  Bight 

Depth 

Cs 

<Cs> 

— 

Cs 

<C5> 

Cs 

<Cs> 

Cs 

<Cs> 

Cs 

<Cs> 

Cs 

<C5> 

(m) 

(m  s‘) 

(ms') 

(m  s~') 

(m  s') 

(ms"') 

(m  s') 

(m  s"') 

(ms") 

(ms-') 

(m  s"') 

(m  s"') 

(m  s"') 

0.00 

1524.9 

0.0 

1537.0 

0.0 

1543.3 

0.0 

1542.8 

0.0 

1525.5 

0.0 

1516.0 

0.0 

10.0 

1524.6 

0.2 

1536.8 

0.1 

1543.0 

0.2 

1542.8 

0.0 

1525.2 

0.2 

1515.8 

0.1 

20.0 

1524.0 

0.5 

1536.4 

0.3 

1542.9 

0.3 

1542,9 

0.0 

1524.9 

0.3 

1515.8 

0.1 

30.0 

1522.8 

1.0 

1535.5 

0.7 

1542.0 

0.6 

1540.3 

0.8 

1524.4 

0.6 

1516.0 

0.1 

50.0 

1519.1 

2.7 

1531.4 

2.4 

1540.4 

1.4 

1537.7 

2.5 

1523.3 

1.2 

1517.1 

0.4 

75.0 

1515.4 

5.0 

1528.1 

4.7 

1538.4 

2.6 

1536.6 

3.9 

1520.4 

2.4 

1512.7 

1.0 

100.0 

1513.1 

6.9 

1525.3 

6.6 

1535.1 

4.0 

1533.9 

5.1 

1517.7 

3.8 

1510.4 

2.4 

125.0 

1511.4 

8.3 

1523.0 

8.2 

1532.9 

5.5 

1529.6 

6.7 

1515.7 

5.2 

1510.1 

3.4 

150.0 

1510.5 

9.5 

1521.2 

9.7 

1530.3 

6.9 

1525.7 

8.7 

1513.7 

6.5 

1508.6 

4.1 

200.0 

1508.0 

11.4 

1519.9 

11.7 

1523.5 

10.2 

1519.7 

12.6 

1509.5 

9.0 

1506.7 

5.5 

Table  lb:  Sound  speeds  and  RMS  sound  speed  deviations  relative  to  the  surface  sound  speed 
as  used  in  the  calculations  described  in  section  3.1.  Data  is  for  the  month  of  July. 


Tasman  Sea 

off  Brisbane 

Coral  Sea 

Exmouth  Gulf 

off  Perth 

Great  Aust.  Bight 

Depth 

Cs 

(m  s"') 

<C5> 

Cs 

<Cg> 

Cs 

<Cs> 

Cs 

<C5> 

Cs 

<Cs> 

Cs 

<C5> 

(m) 

(m  s  ') 

(m  s') 

(m  s-') 

(m  s'') 

(m  s'') 

(m  s'') 

(m  s'') 

(m  s-') 

(ms') 

(m  s'') 

(m  S'') 

0.00 

1511.8 

0.0 

1524.8 

0.0 

1535.8 

0.0 

1534.7 

0.0 

1519.6 

0.0 

1509.7 

0.0 

10. 0 

1511.8 

0.0 

1524.9 

0.1 

1535.7 

0.1 

2534.7 

0.0 

1519.7 

0.1 

1509.5 

0.1 

20.0 

1511.9 

0.0 

1524.9 

0.1 

1535.6 

0.1 

1534.8 

0.0 

1519.8 

0.1 

1509.6 

0.1 

30.0 

1511.9 

0.1 

1525.0 

0.1 

1535.7 

0.1 

1534.8 

0.1 

1519.7 

0.1 

1509.6 

0.1 

50.0 

1511.9 

0.1 

1525.0 

0.2 

1536.3 

1.2 

15334.8 

0.1 

1519.4 

0.1 

1510.3 

0.2 

75.0 

1511.7 

0.1 

1526.6 

0.6 

1536.2 

0.3 

1535.8 

0.4 

1518.7 

0.4 

1510.1 

0.3 

100.0 

1511.4 

0.1 

1526.4 

1.0 

1531.8 

1.1 

1533.2 

0.5 

1517.9 

0.7 

1509.6 

0.3 

125.0 

1510.7 

0.4 

1524.2 

1.0 

1531.1 

2.2 

1530.2 

1.5 

1516.7 

1.2 

1508.6 

0.4 

150.0 

1511.7 

0.4 

1524.1 

0.9 

1527.8 

3.3 

1526.6 

2.9 

1515.2 

1.9 

1506.5 

1.0 

200.0 

1509.5 

0.8 

1522.7 

1.1 

1521.7 

6.3 

1519.8 

6.4 

1510.7 

3.8 

1506.5 

1.8 

In  each  case  the  CR  model  is  compared  to  the  SR  model  based  on  the  surface  sormd 
speed.  The  comparison  is  between  the  predicted  frequencies  corresponding  to  an 
interference  minimum  for  a  series  of  source  ranges  and  depths,  at  several  set  receiver 
depths.  Figure  2  shows  a  typical  LM  fringe  pattern  for  orders  n=0,  1.  The  solid  line 
shows  the  fringes  predicted  using  the  CR  model,  with  the  dashed  line  corresponding 
to  the  SR  model  fringes.  A  single  layer  linear  soimd  speed  profile  (gradient  y  =  -0.1 
ms“Vm  and  Cg  =  1500  ms“^)  was  used.  Clearly  the  CR  model  results  in  a  fringe 
system  of  greater  curvahure  than  the  SR  model,  for  which  the  fringes  are 
asymptotically  linear,  and  at  some  range  (or  time)  the  CR  n=0  fringe  will  overlap  the 
SR  n=l  fringe.  This  motivates  the  use  of  the  fractional  difference 
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rCli  _  fSK 

^  _  Jo  Jo 

f  SR  _  rSR 

J\  Jo 

as  a  measure  of  the  similarity  of  the  two  model  predictions,  where  the  model 
frequencies  are  in  each  case  determined  using  the  same  parameter  set.  For  the 
curvature  of  the  rays  to  have  little  effect  on  the  LM  pattern,  ri=0.  With  r|=0.5,  the  CR 
model  minimum  of  order  n  matches  the  SR  model  maximum  of  order  (n+1),  indicating 
a  significant  deviation  of  the  CR  model  predictions  from  the  SR  model  predictions. 
Hence,  as  a  measure  of  the  dissimilarity  of  the  two  model  predictions,  a  value  ti>0.25 
is  assumed.  It  is  noted  that  T]  describes  the  differences  in  the  shapes  of  the 
frequency/ time  curves  arising  from  the  two  models. 


Figure  2:  Surface  reflections  Lloyd's  mirror  destrucive  interference  fringes  for  orders  n=0  and 
1.  The  solid  line  represents  the  results  of  a  single  layer  xvith  a  sound  speed  gradient  of 
-0.1  ms-'/m;  the  dashed  line,  a  single  layer  with  a  constant  sound  speed  profile.  In  each  case  the 
surface  sound  speed  is  1500  ms-K  The  system  geometry  is:  CPA  range  =  100  m,  speed  = 

10  ms-',  source  depth  =  100  rn  and  receiver  depth  =  100  m. 
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The  deviation  Ti  as  a  function  of  source  depth  and  source  range  for  several  receiver 
depths  is  shown  in  the  contour  plots,  figure  4.  In  each  case,  the  heavy  solid  line  shows 
the  ti=0.25  contour.  Figure  3  shows  the  corresponding  sound  speed  profiles.  Only 
ranges  to  2000  m  and  the  upper  200  m  surface  layer  are  considered.  Over  this  depth, 
the  profiles  exhibit  RMS  sound  speed  deviations  <Cs>,  with  respect  to  the  surface 
sound  speed,  of  ~5— 10  ms”^  in  January  and  —1—6  ms  ^  in  July,  and  maximum  gradients 
of  -0.1-0.2  ms“Vm.  Several  profiles  (notably  in  July)  show  only  small  variations  of 
sound  speed  with  depth. 

It  is  clear  from  figure  4  that  there  is  a  steady  increase  in  the  effect  of  ray  curvature  as 
the  range  is  increased,  as  might  be  expected.  Changes  in  the  source  and  receiver 
depths  result  in  a  less  pronounced  effect,  and  in  several  cases  ray  curvature  is  of  lesser 
importance  with  increasing  depth.  First  consider  the  profiles  with  the  (largest)  steadily 
increasing  gradients  (Exmouth  Gulf,  off  Perth,  Tasman  Sea,  off  Brisbane  &  Coral  Sea 
January  profiles).  All  have  a  similar  form,  and  show  similar  curvature  effects  as 
indicated  by  ri.  Typically,  critical  ranges  (defined  as  the  range  at  which  ti=0.25)  of  600- 
1500  m  are  found,  with  this  range  increasing  as  the  maximum  gradient  decreases.  We 
also  note  that  there  is  a  generally  smooth  evolution  in  the  value  of  Tj  arising  from 
changes  in  the  source  and  receiver  depths.  For  profiles  such  as  those  for  the  Exmouth 
Gulf  and  the  Great  Australian  Bight  in  January,  an  approximate  iso-speed  layer  exists, 
and  expectably  only  when  the  depths  below  this  layer  are  sampled  do  we  see  a  steady 
increase  in  r\  with  range.  A  separate  feature  arising  from  the  geometries  of  the  direct 
and  reflected  rays  is  seen  at  shallow  source  and  receiver  depths.  In  this  case  values  of 
rt  much  larger  than  is  typical  for  the  particular  profile  can  be  obtained  (for  example,  the 
Coral  Sea  January  and  the  Great  Aust.  Bight  profiles).  This  arises  from  the  direct  and 
reflected  rays  having  increasingly  similar  optical  path  lengths,  or  decreasing  time 
differences.  This  will  occur  for  even  small  gradient  profiles,  although  the  matching 
range  will  increase  as  the  magnitude  of  the  gradient  decreases.  At  larger  depths,  this 
feature  will  be  observable  at  larger  ranges. 

The  profiles  can  be  grouped  according  to  whether  or  not  the  sound  speed  has  a 
maximum  at  some  depth  other  than  at  the  surface.  Of  those  that  do  (all  July  & 
Exmouth  Gulf  and  Great  Aust.  Bight  January  profiles),  the  presence  of  local  maxima 
leads  to  shadow  zone  effects  as  well  as  possibly  the  features  already  noted.  This 
shadow  zoning  is  possible  even  in  the  presence  of  small  gradients,  (see  figure  4b  where 
generally  Ti=0).  The  degree  to  which  it  will  be  important  depends  on  the  specific  depth 
and  range  geometry.  In  practice  such  zones  will  not  affect  the  time-frequency 
spectrum  insofar  as  they  will  result  in  a  restricted  range  of  observation,  but  otherwise 
will  not  greatly  affect  the  shape  of  the  fringe  pattern. 

In  as  much  as  these  profiles  represent  a  typical  variation,  evidently  ray  curvature  can 
play  an  important  role  at  close  range.  In  several  cases,  depending  on  the  details  of  the 
sound  speed  profile  and  the  source/receiver  geometry,  the  effects  are  pronoimced  at 
ranges  of  about  600  -  700m. 
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One  measure  of  the  effect  a  particular  profile  will  have  on  the  t]  (Range,  Depth) 
distribution  is  the  "average"  sound  speed  gradient  over  a  given  depth  as  defined  by 

G{D)=<CXD)>ID 

where  <Cs(D)>  is  the  RMS  sound  speed  deviation  at  depth  D.  This  is  motivated  by 
the  reflected  ray  needing  to  sample  the  entire  soimd  speed  profile  between  the  source 
and  the  surface,  and  the  larger  the  variation  in  the  sound  speed  over  this  path,  the 
larger  the  effects  of  ray  curvature  on  the  path  length  (or  transit  time).  From  inspection 
of  the  various  contour  plots,  figirre  4,  and  the  values  G  deduced  from  table  1,  it 
appears  that  a  rough  rule  for  the  adequacy  of  the  SR  model  is  G  <  0.01  ms~^/ m  over  the 
depths  of  interest  (this  corresponds  to  a  nearly  constant  temperature  profile).  This 
value  is  also  suggested  in  calculations  done  using  a  single  layer  of  constant  gradient. 
This  is  for  a  range  to  2000  m,  and  the  value  will  decrease  as  the  range  is  increased.  As 
the  details  of  any  LM  pattern  will  depend  on  the  specific  sound  speed  profile  and  the 
source /receiver  geometry,  such  figures  of  merit  are  of  limited  value,  and  a  direct 
calculation  of  Tj  (or  the  LM  fringe  pattern  itself)  would  be  a  more  productive 
assessment  of  the  adequacy  of  the  SR  model  in  any  particular  case. 


3.2  Data  Comparisons 


Three  sets  of  time-frequency  data,  previously  reported  in  the  AJAAC  bulletin^'^*,  are 
considered.  In  each  case  a  well-developed  Lloyd's  mirror  pattern  is  evident  and, 
owing  to  the  large  water  depths  at  each  location,  can  be  ascribed  to  surface  reflection. 
Further,  source  speeds  are  known  from  independent  analysis  of  the  observed 
narrowband  frequency  doppler  shifts.  It  is  assumed,  and  the  data  support  this,  that 
over  the  acquisition  period  each  source  was  travelling  at  constant  depth  and  velocity. 

There  are  five  variables  of  interest:  the  receiver  depth,  the  source  depth  and  CPA 
range,  the  time  of  closest  approach  and  the  source  speed,  with  the  latter  three  defining 
the  source  range  as  a  function  of  time.  It  can  be  assumed  that  the  receiver  depth  is 
known,  and  the  aim  is  to  replicate  the  observed  LM  broadband  fringe  pattern  and 
thus  deduce  the  source  variables.  As  with  the  SR  model,  the  CR  model  doesn't  in 
general  uniquely  specify  these  variables  and  to  resolve  the  ambiguity  one  of  the  source 
CPA  range,  depth  or  source  speed  must  be  known.  In  these  examples  the  source 
speed  is  specified,  being  known  from  earlier  analysis.  The  remaining  source 
parameters  have  been  determined  using  the  SR  and  CR  models  via  a  weighted  LS 
optimisation  on  a  selected  number  of  LM  fringes.  A  weighting  designed  to  favour  the 
data  nearest  to  the  CPA  time  is  chosen.  This  will  only  be  of  significance  in  those  cases 
where  the  assumed  model  is  inappropriate.  We  consider  each  case  separately. 


13 


DSTO-TR-0064 


3.2.1  Case  A 


These  data^^  clearly  show  the  inadequacies  of  the  SR  model.  Based  on  the  CR  model 
analysis,  ranges  from  ~1200-3400m,  and  corresponding  to  fractional  deviations  of  0.08 
<T1<1.0,  are  sampled  by  the  data.  As  no  experimental  sound  speed  profile  was 
available,  the  appropriate  profile  was  obtained  from  the  TESS^^  database,  and  is  shown 
in  figure  5a.  Based  on  the  receiver  type,  the  receiver  depth  was  assiuned  to  be  305m. 
Table  2  summarises  the  "best  fit"  parameters  for  each  model.  The  corresponding  LM 
fringe  patterns  are  shown  in  figirre  6. 

Clearly,  the  data  show  fringes  with  significant  curvature,  and  that  whilst  the  SR 
model  is  incapable  of  reproducing  this,  the  CR  model  well  matches  the  observed 
fringe  shape.  Furthermore,  the  resulting  "best  fit"  parameters  are  markedly  different 
for  the  two  models.  The  use  of  the  SR  model  in  this  circiunstance  would  have 
important  consequences  for  the  analysis  of  any  narrowband  lines  in  the  time- 
frequency  spectriun.  If  the  imderlying  intrinsic  intensity  profile  of  such  a  line  were 
desired,  then  LM  interference  corrections  and  losses  based  on  the  SR  model  would 
be  inappropriate  and  lead  to  substantial  errors.  If  in  addition  several  such  lines  were 
considered  together,  inconsistent  results  for  the  source  location  would  result  because 
of  the  incorrect  placement  of  interference  minima  at  the  different  frequencies. 


3.2.2  Case  B 


These  data^^  represent  another  example  of  the  inappropriateness  of  the  SR  model. 
Ranges  based  on  the  CR  model  "best  fit"  parameters  cover  ~1000-2400m,  with  0.1<ti< 
1.4.  Again,  the  receiver  depth  is  assumed  to  be  305m.  In  this  example,  a  sound  speed 
profile  recorded  from  a  nearby  location  and  time  was  available^*^  This  is  shown  in 
figure  5b,  whilst  the  "best  fit"  parameters  are  included  in  table  2.  The  LS  fit  results  are 
shown  in  figure  7.  It  is  again  evident  that  the  SR  model  cannot  reproduce  the  data, 
with  the  calculated  fringes  showing  too  little  curvature.  The  CR  model  gives  an 
excellent  representation  of  the  data.  As  before,  the  "best  fit"  parameters  for  each 
model  are  significantly  different,  and  the  comments  above  concerning  narrowband 
lines  are  relevant. 
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Table  2:  Weighted  least  squares  fit  parameters  for  cases  A,  B  &  C.  Parameter  sets  for  the 
curved  (CR)  and  straight  line  (Sl^  ray  models  are  shown.  In  each  case  the  receiver  depth  and 
source  speed  were  fixed  at  the  indicated  values.  The  resulting  curves  are  shown  in  figures  6-10. 


parameter 

Ce^  A 

CaseB 

Cas 

eC 

CR*’ 

SR 

CR*’ 

CR*’ 

SR 

CR*’ 

SR 

CPA  time*’ (s) 

372.5 

372!3 

407.9 

406.3 

405.3 

3573 

357.4 

CPA  range*"’  (m) 

1244.3 

945.2 

958.2 

733.3 

566.8 

823.0 

712.4 

Source  depth  (m) 

198.0 

144.2 

162.2 

120.4 

96.4 

74.7 

62.7 

Source  speed*’  (m  s'*) 

11.4 

11.4 

8.1 

8.1 

8.1 

8.7 

8.7 

Receiver  depth"*’  (m) 

305.0 

305.0 

305.0 

305.0 

305.0 

305.0 

305.0 

a)  time  offset  is  arbitrarily  chosen 

b)  horizontcJ  range 

c)  set  at  value  determined  from  separate  doppler  shift  analysis  of  narrow  band  lines 

d)  depth  unknown,  but  based  on  receiver  type 

1)  fit  based  on  TESS  sound  speed  profile 

2)  fit  based  on  measured  soimd  speed  profile 

Case  A  Case  B  Case  C 


Sound  Speed  (m/s) 

Figure  5:  Sound  speed  profiles  used  in  the  CR  model  calculations  of  section  3.2.  The  heavy  line 
shows  the  profiles  extracted  from  the  TESS  data  base  and  the  light  line  in  (b)  shows  a  profile 
measured  close  to  the  location  of  the  case  B  data  acquisition.  Only  the  surface  350  m  are 
shown. 
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Figure  6:  Surface  reflection  Lloyd's  mirror  interference  minima  contours  for 
The  data  are  described  for  case  A,  section  3.2.1.  The  heavy  lines  show  the  LS  O] 
data,  based  on  the  TESS  profile  of  fig.  5a  (upper  figure)  or  a  constant  sour 
(lower  figure).  The  resulting  LS  parameters  are  given  in  table  2. 
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Figure  9:  Surface  reflection  Lloyd's  minor  interference  minima  contours  for  orders  n=0-8.  The 
data  are  as  described  for  case  C,  section  3.2.3.  The  heavy  lines  show  the  LS  optimal  fit  to  the 
data,  based  on  the  TESS  profile  of  fig.  5c  (upper  figure)  or  a  constant  sound  speed  profile  (lower 
figure).  The  resulting  LS  parameters  are  given  in  table  2. 
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It  is  also  possible  to  test  the  effect  of  using  a  different  sound  speed  profile;  in  this 
case  the  profile  obtained  from  the  TESS^^  database.  Given  that  sufficient  sensitivity 
exists  for  a  discernible  difference  between  the  SR  and  CR  model  results,  it  not 
unreasonable  to  expect  that  the  results  will,  to  a  lesser  degree,  be  sensitive  to  the 
specific  sound  speed  profile.  This  profile  is  also  shown  in  figure  5b  by  the  heavy  solid 
line,  where  it  can  be  seen  that  only  qualitative  agreement  to  the  "true"  profile  is 
obtained.  Lloyd's  mirror  calculations  based  on  this  latter  profile  are  shown  in  figure  8. 
Now  only  qualitative  agreement  with  the  data  can  be  obtained  using  this  soimd  speed 
profile,  and  the  resulting  "best  fit"  parameters  (included  in  table  2)  do  not  agree  with 
the  earlier  values  based  on  the  measured  profile.  Differences  in  the  source  depth  and 
the  CPA  range  of  =25%  are  noted.  The  quality  of  this  fit,  lies  between  that  for  the  SR 
model  and  the  CR  model  based  on  the  measured  profile. 


3.2.3  Case  C 


This  last  example^^  is  interesting  because  the  LM  fringe  pattern  is  not  present  for  the 
larger  source  ranges.  Based  on  the  CR  model  calculations,  ranges  of  —800— 1900m  are 
sampled,  and  we  find  0.06<ri<0.3.  As  for  case  A,  the  sound  speed  profile  is  obtained 
from  the  TESS^^  database,  and  is  shown  in  figure  5c.  The  resulting  LM  fringe  patterns 
based  on  the  two  models  are  shown  in  figure  9,  with  the  "best  fit"  parameters  included 
in  table  2. 

As  suggested  by  the  smaller  values  of  T],  the  SR  and  CR  models  are  equally  capable  of 
reproducing  the  observed  fringe  system  shown  in  figure  9.  We  note,  however,  that  the 
"best  fit"  parameters  derived  using  the  two  models  differ  by  —15%.  Indeed,  as  the 
value  of  Tj  is  not  dependent  on  the  interference  fringe  order,  a  small  value  can  still  lead 
to  a  large  discrepancy  between  the  matches  of  the  two  models  to  the  data  at  higher 
orders  since  any  difference  is  amplified  by  the  order  number.  This  is  shown  in  figure  10 
where  the  CR  and  SR  model  curves  are  calculated  using  the  same  parameter  set  (the 
Case  C  CR  model  parameters  of  table  2).  At  the  CPA  time,  ri=0.06  but  for  the  n=8 
fringe,  a  large  mismatch  between  the  two  models  is  found.  Thus,  although  the  value  of 
Tj  can  give  an  indication  of  whether  the  SR  model  will  be  able  to  match  the  data,  it  will 
not  necessarily  indicate  the  suitability  of  the  deduced  source  location. 


4.  Conclusions 


The  consequences  of  including  realistic  soimd  speed  profiles  on  the  Lloyd  s  mirror 
interference  pattern  have  been  considered,  and  the  need  to  use  such  profiles  has  been 
demonstrated.  Using  ray  tracing  it  is  found  that  ray  curvature  arising  from  a 
piecewise  linear  sound  speed  profile  can  cause  significant  changes  in  the  frequencies  of 
the  LM  interference  extrema  for  a  given  source/receiver  geometry,  and  that  such 
effects  can  occur  within  the  ranges  of  interest  for  passive  sonar  detection.  It  is  also  the 
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being  considered,  together  with  the  source  and  receiver  depths.  Except  where  the 
sound  speed  gradient  is  small  (<0.01  ms"Vm  for  ranges  <2000m)  ray  curvature  is 
potentially  important.  Ranges  of  600  -  700m  have  been  found  to  be  sufficient  to  cause 
marked  deviations  from  the  constant  sound  speed  model.  By  considering  a  fractional 
deviation  frmction,  it  is  possible  to  visualise  the  source  range-depth  regions  that 
potentially  will  be  significantly  affected  by  ray  curvature. 

Three  particular  sets  of  data  for  sources  with  fixed  velocities  and  depths  were 
considered  .  Each  was  chosen  because  clear  broadband  fringe  systems  are  present  for 
times  that  correspond  to  large  source /receiver  ranges,  and  thus  more  readily  elucidate 
any  effects  of  ray  curvature.  It  is  more  usual  that  such  a  broadband  pattern  is  not 
present,  however  as  mentioned  in  the  introduction,  the  effects  of  Lloyd's  mirror 
interference  can  and  often  will  be  present  in  the  (stronger)  narrowband  signals.  These 
data  demonstrate  the  inadequacies  of  the  straight  line  propagation  approach,  whereas 
each  set  can  be  explained  with  the  ray  curvature  model.  It  is  also  found  that  the 
inferred  source  location  is  very  dependent  on  the  choice  of  model.  Thus,  it  is 
important  to  ensure  the  suitability  of  the  applied  model,  for  example,  before  the 
theoretical  LM  pattern  is  used  in  estimating  either  broadband  or  narrowband  sound 
pressure  levels.  A  poor  choice  will  result  in  a  poor  estimation  of  the  intrinsic  levels. 
One  example  is  noted  where  although  a  good  fit  to  the  interference  fringes  was 
obtained  with  the  straight  line  model,  the  resulting  source  location  differs  between  the 
two  models  considered.  The  implication  is  that  care  needs  to  be  exercised  in 
ascertaining  the  applicability  of  a  constant  sound  speed  profile. 

For  a  further  example,  it  was  possible  to  compare  the  results  of  calculations  using 
two  different  sound  speed  profiles  with  the  same  acoustic  data.  The  results  indicate 
that  the  quality  of  any  fit  and  the  thus  derived  source  parameters  will  be  sensitive  to 
the  choice  of  profile.  It  is  therefore  important  that  a  well  described  true  profile  (ie.  one 
measured  temporally  and  spatially  close  to  the  true  location)  be  used.  Given  the 
number  of  potential  unknowns,  this  also  suggests  that  the  receiver  depth  be  correctly 
known,  as  the  specific  source/ receiver  geometry  is  important  in  determining  the  effect 
of  ray  curvature  for  a  given  sound  speed  profile.  Information  that  can  be  obtained 
from  an  independent  analysis  of  narrowband  lines  {viz.  source  speed,  approximate 
CPA  range)  is  also  of  use,  given  that  a  unique  solution  is  not  guaranteed  from  the 
Lloyd's  mirror  analysis  alone.  With  this,  it  should  be  possible  to  correctly  utilise  the 
information  inherently  present  in  the  time-frequency  spectrum. 
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A  knowledge  of  the  interference  between  direct,  surface  and  bottom  refelected  rays  emanating  from  a  moving 
underwater  source  can  be  used  to  determine  the  source-detector  relative  positions.  Using  simple  ray  tracing,  the 
frequency-time  interference  patterns  for  a  broadband  frequency  source  arising  when  the  propagating  medium 
has  a  piecewise  linear  souiad  speed  profile  are  considered,  and  are  compared  with  those  for  a  constant  sound 
speed.  It  is  found  that  even  small  deviations  from  a  constant  sound  speed  can  lead  to  a  noticeably  different 
interference  fringe  placement  and  shape,  and  consequently  to  the  inferred  source-detector  geometry.  Further, 
comparisons  of  model  calculations  to  experimentally  observed  fringe  systems  demonstrate  the  inadequacy  of  the 
constant  sound  speed  model,  whereas  the  use  of  the  simple  segmented  profile  is  capable  of  replicating  the 
observed  patterns. 
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